Improvement of the matching of the exact solution and variational approaches in an 

interacting two-fermion system 
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A more reasonable trial ground state wave function is constructed for the relative motion of 
an interacting two-fermion system in a ID harmonic potential. At the boundaries both the wave 
function and its first derivative are continuous and the quasi-momentum is determined by a more 
practical constraint condition which associates two variational parameters. The upper bound of the 
ground state energy is obtained by applying the variational principle to the expectation value of the 
Hamiltonian of relative motion on the trial wave function. The resulted energy and wave function 
show better agreement with the analytical solution than the original proposal. 
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I. INTRODUCTION 

Tremendous progress in cooling, trapping and manipu- 
lating ultracold samples makes it possible to experimen- 
tally study few-body phenomena in trapped atomic and 
molecular systems with precise control The study 
of trapped few- atom systems provides a bridge between 
small and large systems and the treatment of such meso- 
scopic systems is complementary to the interpretation 
derived within effective many-body frameworks. Systems 
made of just two interacting atoms are especially relevant 
due to its role as building blocks of many-body strongly 
correlated states. Two-body system with tunable inter- 
actions has been realized in the Heidelberg experiment 
using two fermionic ^Li atoms in the ground state of 
a potential created by an optical dipole trap. In this ex- 
periment, the interaction energy of two distinguishable 
fermions as a function of the interaction strength has 
been measured showing the experimental capability to 
simulate strongly correlated few-body quantum systems. 
An exceptional property of the two-body interacting sys- 
tem is that there exists analytical solution for arbitrary 
values of the interparticle interaction [1, 13) ■ A theory for 
the tunneling of one atom out of a trap containing two 
interacting cold atoms is developed by introducing the 
quasi-particle wave function Q. Unfortunately, no the- 
oretical solution is available for more complex systems 
even if one more atom is included. 

It is thus of interest to seek for some approximation 
methods which may potentially be generalized to more 
than two particles. Ref. @ provides such an example, 
which calculated the ground-state energy of the relative 
motion of a system of two fermions with spin up and spin 
down interacting via a delta-function potential in a one- 
dimensional (ID) harmonic trap. The authors divided 
the relative coordinates space into three regions and es- 
tablished a trial ground state wave function according 
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to the Bethe-ansatz solution for a system with delta- 
function interaction and the ground state wave function 
of harmonic potential. They directly applied the periodic 
boundary conditions to obtain the Bethe ansatz equation, 
which was taken as the constraint condition for quasi- 
momentum. By analyzing solution of the Bethe ansatz 
equation, we know that the possible values for the quasi- 
momenta k depend on the value of c 0] . For the repulsive 
case (c > 0), only real values of k are ground-state so- 
lutions of the Bethe ansatz equation. For the attractive 
case(c < 0), only the purely imaginary numbers k satisfy 
for the ground state. For both cases they found the nu- 
merical minimization of expectation value of the ground 
state energy by introducing two variational parameters. 
This work provides a starting point for the investigation 
of more complex few-body systems where no exact theo- 
retical solution is available. 

In this paper we improve the trial ground state wave 
function proposed in [6| from two aspects. Firstly not 
only wave function but also its derivative are continu- 
ous at the boundaries, which avoids introducing the dis- 
continuity of the first derivative of the wave function for 
regular potential. Secondly instead of using Bethe ansatz 
equation as the constraint condition, we find a more prac- 
tical constraint condition for the quasi-momentum which 
incorporates the two variational parameters and the con- 
tact interaction strength. There, the solutions of the con- 
straint condition no longer depend on the sign of c. 

The paper is organized as follows. In Sec. II, we in- 
troduce the modified trial ground state wave function for 
the system of two interacting fermionic atoms. At the 
boundaries both the wave function and its first deriva- 
tive are continuous and we obtain a similar constraint 
condition for the two variational parameters. Then the 
variational principle is applied to the expectation value 
of the Hamiltonian of relative motion on the conjectured 
wave function in Sec. HI. In Sec. IV Our results are com- 
pared with the analytical result obtained in 's', and the 
approximate result obtained in Q. In addition the den- 
sity distributions for various interaction constants also 
show good agreement with the analytical solution. We 
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conclude our results in Sec. V. 



II. TRIAL WAVE FUNCTION 

We shall consider the same system of two interacting 
fermionic atoms with mass m loaded in a harmonic trap 
of frequency ui as in Ref. Q. The Hamiltonian can be 
formulated as 



ft2 d' 



where Va {xi,X2) is the trapping potential for atoms lo- 
cated at xi and X2 

Va {xi, X2) = —muj^x\ + —mLo^x^. (2) 

The interaction potential V/ is assumed to be s-wave con- 
tact potential represented by a delta-function 



Vi {xi,X2) = 2c5 (x2 - a;i) 



(3) 



where the interaction strength c can be tuned via mag- 
netic Feshbach resonance from —00 to +00. 

Introducing the center-of-mass and relative coordi- 
nates X — xi~X2,X — {xi+X2)/2 allows us to decompose 
the Hamiltonian ([T]) into the center-of-mass part 



HcM — 
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2M dX^ 
and relative motion part 



(4) 



(5) 



The Hamiltonian of the center-of-mass is nothing but a 
simple harmonic oscillator with total mass M = 2m, 
while the relative part describes an oscillator with re- 
duced mass fj, — m/2 and a delta potential well/barrier 
in the trap center. The latter has been exactly solved in 
Ref. 0. 

In Ref. the authors built a variational trial wave 
function iIj{x, {a, L}) based on the Bethe ansatz solu- 
tion of two particles and it was shown that applying the 
same principle for a large number of fermions is possible. 
The parameters a and L controls the decay of the trial 
function outside the trap and the region where the decay 
occurs. Minimizing the expectation value of the Hamilto- 
nian Hrei on the conjectured wave function ip{x, {a, L}) 
by setting the derivatives with respect to these two pa- 
rameters zero gives the upper bound of the ground state 
energy of the system 



Egs < 



(6) 



The trial wave function ip is continuous in Ref. [6| , how- 
ever, the first derivatives of the wave function develop an 
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FIG. 1: The potential and the trial ground state wave function 
4' (a;) in the relative coordinates system of interacting two- 
fermion system. Both the wave functions and its derivatives 
are continuous at the boundary x = ±L. 



unphysical discontinuity at a; = ±L. We know from basic 
quantum mechanics that in matching the wave functions 
in different regions the wave function -0 is always contin- 
uous and the first derivative dtjj/dx is continuous except 
at points where the potential is infinite, i.e. S poten- 
tial. In our problem there do not exist such S potential 
points at a; = ±L hence discontinuity of the first deriva- 
tive there should be excluded. As we shall see below a 
simple improvement of the trial wave function will help 
us get rid of this dilemma and the resulted ground state 
energy matches the analytical result even better. 

In the following we introduce the variational ansatz 
and the modification of the trial ground state wave func- 
tion. As shown in Fig. [1] the wave function is continuous 
and it's derivative is also continuous, which is much closer 
to the real wave function. 

How do we build the trial ground state wave function 
for the two body problem of Hamiltonian i?rei? The 
space is divided into three regions by the parameter L. 
The wave functions in the regions I and HI are chosen 
as eigenstates of ID harmonic oscillator center at a; = 0, 
instead of at x = ±L as done in Ref. Q. "0/ 
have the form of the Gaussian with the parameter a. 
In the regions II, i.e. in the vicinity of the S potential 
center, the contact interaction term plays the dominant 
role and the harmonic potential can be neglected. The 
wave function takes the form of two distinct fermions 
with a contact interaction. Formally we assumes the trial 
ground state wave function as the following configuration 



ip = < 



' ipj = Ae"""'-^ -00 < X < -L 

^Pjj = (Be-''^^ + Ce'''^) e (x) 

+ (De''^^ + Ee-''''') 8 (-x) -L < x < L 
i^iii = Fe'"^^ L < X < 00 

(7) 
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where Q is the Heaviside step function 



(8) 



The remaining work is to match the wave functions 
in different regions through the boundary conditions. In 
region II, the continuity of the wave function ipjj at a; ~ 
gives the fohowing coefficient relation 



B + C ^ D + E. 



(9) 



By integrating the eigenequation Hipn = Eijjjj from the 
negative infinitesimal to the positive infinitesimal, we do 
find a jump of the first derivative of wave function at 
x = 



ik{-B + C ~ D + E) 



X 2c{B + C). 



(10) 



Because the potential is axially symmetric, the ground 
state wave function must be spatially symmetric and has 
zero nodes, so 



B = D,C = E,A = F. 
According to above three equations, we get 

ik 



B = 



ik ■ 



2tic 



2/ic 



(11) 



(12) 



Kix = —Lovx — L^ the wave function is continuous, 
which means 



Ae 



C 



■ikL 



(13) 



and the first derivative of ■(/' is continuous too, we get 



2aLe-"^ A = Cik 



ifc + -p- 



JkL 



(14) 



The combination of (|13p and (|T4|) leads us to a constraint 
condition 

'^ck + i2La 



i2kL 



k + i-. 



i^r k - 



i2La 



(15) 



This equation plays similar role as that of Bethe ansatz 
equation. It is, however, irrelevant with the boundary 
conditions (periodic or open) for the two-fermion sys- 
tem which should always be chosen before one deals with 
the exactly solvable models. In this way we avoid ar- 
tificially imposing the periodic boundary condition onto 
the system. And we find no obstacles to generalize this 
constraint condition to N fermions following a similar 
scheme as in the Appendix of Ref. Q . 

On the basis of above restrictions to the coefficients, we 
get the following trial ground state wave function subject 
to normalization 



tfjjj — cos{9 — kx)Q (x) 
+ cos{0 + kx)Q (-x) 
III = ^//(i)e-«-'e«^' 



— oo < a; < — i 

-L < X < L 
L < X < oo 



(16) 
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FIG. 2; (Color online)The ground state energies of the rel- 
ative motion for two distinct fermions interacting via s-wave 
pseudopotential confined in an axially symmetric harmonic 
trap with angular frequency uj. The black dotted lines indi- 
cate energy levels of analytical results, (a) The blue curves 
represent the energy by using Bethe ansatz equation as con- 
straint condition, (b) The red curves represent the energy by 
using constraint condition of the continuity of wave function 
derivative. 



where 9 satisfies the following relation 



9 = arctan 



kti^' 



(17) 



This trial ground state wave function, composed of pure 
elementary functions, offers convenient analysis due to 
the continuity of both its first derivative and itself We 
schematically show the configuration of the potential and 
the wave functions in Figure [TJ which join each other 
smoothly at the boundaries for all three regions. 



III. VARIATIONAL PRINCIPLE 

We choose a and L as our variational parameters and 
aim to find the least possible value of the ground state 
energy on the trial wave function for different interaction 
strength c subject to the constraint condition ([T5|) . Con- 
trary to the original proposal (eq. (3.2) in [6|), in our 
scheme, both variational parameters a and L enter into 
the constraint condition, which determines the quasi mo- 
mentum k. To apply the variational principle, we need to 
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compute the normalization factor of the wave function, Moreover, the expectation value of H^ei on the ground 
which yields state pB]) is calculated as 

l + cos(2fcL-20) ri /tttx^ 

^"^'"^^ " V8^ (1 - erf(V2ai)) 

+ ^(sin(2fcL-26i)+sin20 + 2fcL) (18) 
2k 

where the error function erf(a;) is defined as 

erf(a;) ^ —= e^^'dt. (19) 
A Jo 

I 



I Hrei = [1 + cos(2fcL - 261)] {Aa^h^ + n^uj^^ 



TT l-erf(2aL)) , L \ 



2a 16AiQ;e-2"i' 8fia J 



[sin (2kL - 29) + sin (261)1 + '■ — + + c (1 + cos 261) 

4/i ^ ^ ' ^ " 2iJ. 6 

2 

[{2k'^L^ - l) sin {2kL - 261) + 2fcL cos {2kL - 20) - sin 26l] 



(20) 



It is interesting to consider first the limiting case of 
L — !■ 0,c — )■ 0, that is, a simple harmonic oscillator. Ac- 
cordingly the upper bound for the ground state energy 
dni) reduces to 



lim 



rel 



2^i 



IIUJ 

~8a 



(21) 



The minimum value of the energy can be found by setting 
the first derivative with respect to a to zero. We find the 
extreme value of a, which indicates the decay of the trial 
function, to be a* = ^ui/2h, such that 



hm , ,1 ,. 

L,c-^0 



1 



(22) 



In this case we recover precisely the ground state energy 
of the simple harmonic oscillator. 



IV. COMPARISON 

For the model system the constraint condition 
the trial wave function reduces to 



k'^h? tan {kL) — 2cfik 
2h'^Lk + 4^icL tan (kL)' 



51) for 



(23) 



Putting this back into ([6]) , we sweep over all values of L 
and k for each assigned c and minimize the expression ^ 
numerically. We call the extreme parameters L* and k* 
when the minimization point has been reached. In this 
way we determine the ground state energy of the two- 
fermion system as a function of the coupling c via the 
variational principle, for both repulsive and attractive 



interaction. The result is depicted in Fig. [5] in terms of 
the physical variables e and aiD , which are related to the 
energy Egs smd interaction strength c through 



E, 



GS 



aiD = 




(24) 



and the value of e is determined by the follow implicit 
equation 



2aiD 



r(- 



r(-f) 



(25) 



where T{x) is the complete gamma function. We notice 
that by combining the variational principle and continu- 
ity of the wave function derivatives a better agreement 
is achieved between this analytical result (black dotted 
line in Fig. [5]) and our result (red solid line in Fig. [^b) 
than that obtained with the Bethe ansatz equation as 
constraint condition (blue solid line in Fig. [5^). In the 
repulsive case the modified variational curve lies exactly 
on the top of analytical one and the effective matching 
range of the approximation in the attractive case is also 
much wider than the original proposal. 

Inserting the parameters L* and k* back into the equa- 
tion (1231) I we get the corresponding a* . The trial wave 
function pB)) is determined by these parameters L* , a* 
and k* for each value of c. The probability densities of 
the trial ground state wave function in the relative coor- 
dinate are illustrated in Fig. [3] over the whole range of 
c. It reveals that for repulsively interacting system our 
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FIG. 3: (Color online)The normalized probability density in 
the relative coordinates of two interacting fer mions. Here 
the inter-particle interactions are in units of ^/ ^/h?u) and 
the ID scattering lengths are in units of y/hfjujo-. (a) c — 
0.05, aiD = -10; (b) c = 0.25, am = -2; (c) c = 20, aio = 
-0.025; (d) c = -0.5, aio = 1; (e) c = -0.1, aio = 5; (f) 
c = —0.05, aiD = 10. The red solid lines show the analytical 
solutions of the probability density while the black dotted 
lines show the variational results with modified trial wave 
functions. 

variational wave functions near perfectly match the ana- 
lytical ground state wave functions of the relative motion 
Hamiltonian Hrei (see Fig. [3^-|3J:), which is 

ipo {x) = AoU (^-|, 0.5, 2ax^^ exp {-^iujx^/2h) . (26) 

The extension to weakly attractive case is possible (Fig 
and[3f), but not too far. If the interaction is negative. 



the ground state is McGuires cluster state [1, Q , which 
is a bound state and would decay quickly via molecular 
channels. We see that whereas our variational theory 
yields surprisingly accurate energies, the wave functions 
are notoriously poor for finite attraction (Fig. [3ji)- 



V. CONCLUSION 

In conclusion, we have found a much better trail wave 
function that approximates the ground state energy of 
two interaction fermions in a ID harmonic trap via vari- 
ational principle with the constraint conditions of the 
continuity of both wave function and its derivative. This 
provides a modified scheme which enables us to get rid 
of the idealized periodic boundary condition in the cen- 
tral region —L<x<L. Instead we connect the wave 
functions and derivatives at x — ±L by simple quantum 
mechanics rules. A constraint condition is obtained to 
determine the quasi momentum, under which the varia- 
tional principle is applied to minimize the ground state 
energy. This offers a new approach to solve the few-body 
problem and the better agreement between our results 
and the analytical solution shows that it is convincing to 
extend our method to many-body systems, which arouse 
experimental interest [Toj . 
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